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Abstract 

In this work we consider fermionic zero modes in the external scalar and electromagnetic field 
forming the vortex on a sphere. We find the correspondence between the equations for the fermions 
in different dimensions, find their explicit expressions through the vortex fields in case of massless 
fermions, asymptotics near the poles in case of massive fermions and check the number of the 
solutions by proving index theorem for the fermions on a sphere. As a part of deriving the index, 
we write a detailed calculation of the Green function of the Heat equation. 
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I. INTRODUCTION 



Recently, considerable interest has been attracted to localized fermionic states in topo- 
logically non-trivial external fields. Its revival is stimulated by developments in two very 
different areas of physics, high-energy particle theory and condensed-matter physics (theory 
and experiment). 

The Standard Model of particle physics, though extremely successful in description of 
most of the experimental results, cannot explain the apparent hierarchy between two funda- 
mental scales, the electroweak symmetry-breaking scale v ~200 GeV and the gravitational 
(Planck) scale Mp\ ~ 10 19 GeV. It has been suggested [if] and discussed in thousands of 
research papers (for reviews, see e.g. Ref. j^, 3|) that the gravitational scale may be much 
lower if there exist additional space dimensions beyond the three observed ones. Their 
non-observation may be related either to compactification of extra dimensions down to a 
very small size {4, 5] or to localization of visible particles on the observed four- dimensional 



Minkowski subspace 



, |7| . The latter may naturally happen due to interaction of a particle 
with the external field with non-trivial topological properties (e.g. {{j, 9]). 

Quite recently, it has been realized that a similar situation may occur in a laboratory, 
where non-trivial topology of the effective bosonic external field results in the localization 
of fermionic states at particular points, lines or surfaces in certain materials. The motion of 
fermions is governed in these cases by the very same equations as in the particle-physics mod- 
els and in many cases, the localized states have been observed in a laboratory. Well-known 
examples include the fractional Hall charges [lO], topological insulators [1 if ] . topological su- 
perconductors 12] and graphene [3]. The topological nature of the localized states makes 
them in principle suitable for the purposes of quantum computing jl4|. Reviews of some 
topics in this rapidly developing field may be found in Ref. [15 |. 

In this paper, we will study localized fermionic states in the external field of abelian 
Abrikosov-Nielsen-Olesen vortex [if]. \y\. 

The purpose of the present work is to explore various possible fermion- vortex interactions 
for various number of space dimensions. In (3+1) dimensions, it has been shown that for 
the vortex winding number n, there exist n localized zero modes of fermions with either 
Majorana-like {9] or Dirac-like 18] equations. These zero modes describe (l+l)-dimensional 
massless fermions which move freely in the core of the vortex line dubbed "superconduct- 



ing strin g" ; the existence of these modes is related to the vortex topology by the index 
theorem 19]. This feature has been used for extra-dimensional model building jl[ Q] in 
(5+1) dimensions, where the vortex fields depend on two extra spatial coordinates while 
the core corresponds to our three-dimensional space. Two (related to each other) features 
of the localized modes are important for construction of realistic particle models: chirality 
and masslessness. The former provides for correct quantum numbers of the standard-model 
fermions while the latter opens up a possibility to generate particle masses at a correct scale 
by means of an additional interaction with the Higgs scalar field treated perturbatively. 
Moreover, if three fermionic generations are associated with three zero modes of an n = 3 
vortex, then an elegant mechanism to explain the mass hierarchy among the families of 



charged fermions may operate 
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2l| . It has been shown 2CH22I] that within this approach, 
not only the hierarchical pattern of masses and mixing may be reproduced, but their partic- 
ular experimental values are obtained with a surprisingly small number of free parameters. 
Finally, the very same model may automatically explain [23j the observed hierarchical pat- 



tern of neutrino masses and leptonic mixings whch is very different (e.g. Ref. 24j ) from the 
one observed among charged particles. 

Though various compactification schemes for the extra dimensions may be considered, 



a practical working example [2l| for the quoted phenomenological results was the M 4 x S2 
manifold, where M4 is the usual (3+l)-dimensional Minkowski space and S2 is the extra- 
dimensional sphere. Though the vortex on the sphere has a somewhat different topology 
from the flat space (for instance, nontrivial topological configurations may exist for purely 
gauge fields), it has been shown |2l| that the basic prperties of the flat-space solution remain 
intact. However, to the best of our knowledge, the index theorem for the spherical case has 
not been proven. 

Returning to condensed-matter examples, we point out that the vortex-on-the-sphere 
equations are relevant for description of fermionic states on a fullerene molecule, where the 



25). 



effective vortex-like external field is induced by disclinations in the hexagon vertices 
Possible applications of the fullerene-localized fermions are yet to be explored. 

Motivated by these important applications, we consider here the fermions with Majorana- 
like interaction with the vortex field. The fermions are assumed to live on a M x S2 manifold, 
the vortex field depending on the two S2 coordinates. We will relate Majorana and Dirac 
cases for arbitrary D to each other (and ultimately to the cases of lower dimensions), study 



the properties of the fermionic zero modes and derive the index theorem. 

The rest of the paper is organized as follows. In SecJTT] we write the vortex-fermion in- 
teraction, study the fermion equations of motion and perform the reduction of the cases of 
different dimensionality to each other. We find the structure of the fermionic zero modes by 
a direct analysis of these equations made in|B] In Sec JIIII we derive the corresponding indices 
and prove particular index theorems which of course are in agreement with the explicit results 
of SecJTT] We summarize our conclusions, emphasise properties of two phenomenologically 
most interesting cases (M4 x S2 - extra-dimensional extension of the particle-physics Stan- 
dard Model ) and discuss possible applications of other cases in the section [TV] Notations 
and technical details are collected in the Appendices. 



II. STRUCTURE OF ZERO MODES 

In this paper we consider fermions on multi-dimensional manifold S 2 x M, where M is 
Minkovski space. If the dimensionality of the manifold is even, we designate it by D. If 
the dimensionality is odd, it is convenient to designate it by D + 1. The sphere contains an 
Abelian gauge and a scalar field which form Abrikosov-Nielsen-Olesen vortex. Its structure 



is described in detail in ref. 



211 ] . We briefly recall it here. The Lagrangian for these fields is: 



L = VG~ s (-\F ab F ab + (D a (j)) + D a (j) - ^(0 2 - v 2 ) 2 ), (1) 

where D a (f> = d a <t> — ieA a (p and F a b = d a Ab — dbA a (all kinds of subscripts are described in 
appendix [A]). To describe the solution we have to introduce two patches (this discourse was 
proposed in ). The first of them covers the sphere with the South Pole excluded. In 
this patch one may use the following Ansatz for the solution: 

= /(*)e* 
A v = A(0), 

A e = 0. (2) 

The other patch covers the sphere with the North Pole excluded. In the overlapping region 
its field values are connected to the values on the first patch by the gauge transformation 
with the gauge function equal to (—<f) , 

A — y A—-, — ► 0e~T (3) 
e 



It can be shown that f{6) and A(9) have the following asymptotics: 



6^0: f~C 0F 9, A~C 0A 6 2 , 

e^TT: f~C vF , (4) 

e 

where Cof, Co a, C^f are constants. 

Now let us consider fermions interacting with the vortex. To construct the Lagrangian 
we introduce velbain which is: 

p a — — A a 
~~ R ' 

e a = _ S a 

* Rsm6 *" 

e« = 5« (5) 

Then we define the spin connection R°^ = |e^V b^ A,iS ■ It can be shown that R^ 2 = — -R 2 . 1 = 
and all other components are equal to zero. The Lagrangian for the fermions is defined 

as 

L = V^G (iNT a {D) eiD A N + {U 2k N + T\ D) C {D) N* + /i.e.)) , (6) 

where 

R 0J 

D A = d A + -f-T[ D) Yj D) - ie(k + kT? d) )A a , (7) 

g can be chosen a real coupling constant without loss of generality (if g had a complex 
phase, we could include it into N ), k and k are integer or half-integer charges. So in case 
we consider the Lagrangian is: 



L = v / -G {iNV; D) d,N+ zN^(de + + 

r 2 

iR RM {dip " ieMk + ^ <P$ N + 2 ( $2fciV+r (^)^)iV* + h.c.)) (8) 

Let us note that since N is a column of anticommuting variables, the term responsible 
for the interaction with the scalar field in flS]) is non-zero only if T^ D ^C^ is antisymmetric. 
If we take into account (1A13j) we see that in even number of dimensions a few cases are 



possible: 



• D/2 mod 4 = 0, no term with majorana mass 

• D/2 mod 4 = 1, (a) C {D) = (7$ 

• D/2 mod 4 = 2, (b) C = C^, (c) C = C@ } 

• D/2 mod 4 = 3, (d) C = Cf^ 

Also if we consider a gauge transformation of the term |$ 2fc N + T°^C(e>)N* , we can see 
that the term containing k is allowed by the gauge symmetry only if the anticommutator 
{r^)) r^CfD)} = 0, that is in cases (a),(d). 

If the number of dimensions is odd and equal to D + 1, to obtain Lorentz-invariant 

(2) 

couplings we cannot assume C = C)jL or write the term containing k. So we are left with 
the following cases: 

• D/2 mod 4 = or 3, No term with Majorana mass 

• D/2 mod 4 = 1, C (D) = 

• D/2 mod 4 = 2, C [D) = cfy 

Thereby the equations for odd number of dimensions can be considered as a special case of 
the equations for even number of dimensions except for D = 2. 
The Lagrangian results in the following equations of motion: 

r 2 

'IT^ - ieA v(k + «f f D) ))N + g<S> 2k C (D) N* = (9) 
ii sin u 

N can always be decomposed into a set of plane waves propagating in M. Separating positive 
and negative frequency, we get the following expression: 

N{x A ) = J d A k{N + {k^x a )e ik ^ + N_{k^x a )e- tk ^). (10) 

After substituting this decomposition into eq.([9]), multiplying it by T and separating the 
terms with different exponents we have the following equations: 

(-k - T {D) T\ D) ki)N + + D kin N+ + g& k T\ D) C {D) N*_ = 0, 
(k + r° {D) T\ D) ki)N- + D km N. + g® 2k T\ D) C {D) N* + = 0, (11) 



where 

T {D) L {D) /a COttf . L (D) L (D) 



D tln = i T -^Si( Ss + ^) + i^f^a, - + «f f D) )). (12) 

it 2 Rsmv 

We introduce a bispinor 



A = ( N+ | (13) 

AT* 

and rewrite the equations of motion as 

A;°A> = i^W + DN, (14) 



where 



r° r* n 



and 



d=[ Dkm ^Vfa] (15) 

V -^(^*) 2fc (r? D) C(D))* -D* kin J 

One can note the following anticommutation property of A 1 : 

= 0. (16) 

By making use of the last relation, we conclude that if iV is an eigenvector of D with an 
eigenvalue m: 

DN = mN, 

then Na = where \k\ = k l k l , is also an eigenvector of D but with opposite eigenvalue, 

DN A = -mN A . 



Also one can see that if N satisfies Eq. ffl4l) . Na also satisfy this equation. These relations 
can be written in the following form: 



This system has non-trivial solutions when (A; ) 2 = (k 1 ) 2 + m 2 . In this paper we will 
explore the zero modes, that are eigenvectors with m = 0. 

So let us consider the equation DN = which can also be written as 



Dun -^^ + g& k T\ D) c {D) ( J'V" J =°- ( 18 ) 

We see that these two equations have similar structure. Further we will show that they 
have a finite number of fundamental solutions . If we find them, we can write a general 
solution as: 

iV + + iV_ 



^-^ = £/W ( °> (19) 
i 

or 

N+ = (a l +ip l )N®, (20) 

iV_ = (a,-iA)^ ( °. (21) 

Here a; and $ are arbitrary anticommutative coefficients. They are real because the equa- 
tions are linear only in terms of multiplying by real numbers. Also let us note that since 
oci,/3i are anticommutative are commutative. From the last formulae we conclude that 
iV + have the same structure as iV*. 

Now let us consider the fundamental solutions. We will reduce the higher- dimensional 
equations to the case of two dimensions. The asymptotics of the solutions are discussed in 
the appendix [B] The equation for the fundamental solutions is: 

g m 2k c [D) (N^y = o (22) 

We note that inside the vortex N"> obeys the effective Dirac equation in M, 

iT^N = 0. 

A 2 £) / 2 -component spinor can be represented CIS db column consisting of two (2 D ' 2 1 )- 
component spinors Ni,N 2 : 

N® =1^1. (23) 
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s 



Let us substitute this expression for N"' into the eq. (T2z]) . The result is different for each 
case so let us consider them separately. 

In case (b) Eqs. (l22|) decouple into two equivalent systems: 



2 



iT\ D _ 2) {d e + ^V 1)2 + i-^(d v - ieA v k)N 1)2 + g m 2k C^_ 2) Nl 2 = 0. (24) 
In case (c), the equations decouple into the systems 

1 r 2 

iT\ D _ 2) (d e + -^-)iVi + - ieApQN! + igR^C^Nl = 0, 

t o r 2 

iT\ D _ 2) {d e + —)N 2 + i-^fi 9 * ~ ^A v k)N 2 - i g m 2k C^_ 2) N; = 0. (25) 
By the phase transformation 

JVi ->• e ln/4 N u N 2 -> e- iw/4 N 2 (26) 

we reduce the equations to the same form as in case (b). 

We can see that Eqs. ( I24j) have the same form as the equations of motion ( 122]) for the 
fermions in (D — 1) or (D — 2)-dimensional space. So both cases (b) and (c) are reduced to 
the case (a). 

Now let us consider D-dimensional fermions in case (a). After the substitution ( 123|) the 
equations become: 

to r 2 

^m-Ade + -=-)JVi + - ^A v (k + k))^ - i 9 m 2k C%_ 2) N* = 0, 

z sin o v 1 

to r 2 

tt\D-2)(de + ^V)iV 2 + i-¥=£-(d v - ieA v (k - k))N 2 + igR$ 2k C%_ 2) N* = 0. (27) 
z sin o K ' 

Let us multiply the second equation by Cp_ 2 ^ and complex conjugate it. Then we introduce 

a new variable 



N 2 = iC$_ 2) N* (28) 



and obtain the following equations: 



cot0,„ r 



2 

zlV 2) (0 fl + ^)N 1 + f-^-(d v - ieA^k + k))N! - g m 2k N 2 = 
v ; z sin v 



cote,- r 



2 

zl%_ 2) (<% + ^)iV 2 + z 1 ^^ + ieAp{k - /?))iV 2 - gR($*) 2k Ni = 0. (29) 
z sin u 

We now introduce a new spinor \1/ defined as 

''' I ^ I (30) 



2 



It describes zero modes of Dirac-like interacting fermions with the Lagrangian 

g^ {1 ~ T ^K - m ^ 1 + T f D ) ) ^ ) (3| , 



and the equations for zero modes 



cote,, r 



2 

L sin f 

giZS 2 ^ 1 "^ * + gRm 2k il+ ^ D)) y = 0. (32) 

In case (d) we again make the change (I28p . Then we repeat the steps made above, make 
the same change (13"U1) and obtain the equations: 



pup* + + - + fcrg,)))* - 

V)^ _ gRm ^llhBl^ = o. (33) 
2 2 



Now let us reduce the system (1331) to the case .D = 4. If we represent \I/ as a set of iV 

_D — 2 

two component spinors . . . ^jy where iV is equal to 2 _ 2 _ : 

and use ( 1A4I) . then the system (1331) decouples into the following set of equations: 

2 

( 2 ) ^ „•„ 4 n„ , „DA2fc, 



^ r (2)(^ + ^H^n + - ieA v (k + - ^$ 2fc *Ar_ n+1 = 



2 7 sine 



cote,, . r 



2 

«rf 2) (^ + ^)^v_ n+1 + - ze^(K - fc))*iv-n+i - gR($*) 2k V n = 0, (34) 

w 2 sine 

n = 1..JV/2. 



One can see that the spinor 




satisfies Eq.([33]) for D = 4. 
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In the same way we can make a reduction from four to two dimensions. We write ty n and 
^N+i-n in the component form: 



* t 




N+l-n 



and after using flA3j) rewrite Eqs. fl34l) as 



(35) 



ir; 2) («9 e + 



COt0 N 




r\ 2) {d v -ie{K + kT\ 2) )A^ 



m9 ® 



r° i -+- r° 

£>$2fc + 1 + 1 (2) 




e 1 



VN+l-n 



0, (36) 



*rf 2) (^ + 




+ ^ r (2)(^-^( K " A;r ( 2 ))^ 



sin 



gR 



1 + 1 (2) $2fc+ jW ^afc 




0. (37) 



One can see that the first of these equations is equivalent to ( 1331) in two-dimensional case, 
and the second one becomes equivalent after the change k — > — k, $ — > $*. In the appendix 
|B]the component form of these equations is given. 

Thus we showed that the D dimensional spinor satisfying the Dirac equation for zero 
modes ( 133]) can be expressed in terms of a 2-component spinor satisfying the Dirac equation. 
The structure of fundamental solution is the following: 




(38) 



where 



I o N 



V / 



^2 



N+l-n 



* 7 




N+l-n 



V N+l-n 
^N+l-n 



n 



1,..., N/2. 



(39) 
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For g 7^ the only non-zero components are CniVN+i-n- ^ 9 = 0; non-zero components are 
linearly independent but their structure is more complicated and it is given in the table H 
Explicit expressions and asymptotics of the solutions are found in the appendix [0 The 
spinor has the following structure: 



case (a), 



*»=\ J; h (40) 



case (b), 



N H) = | <} | j Nil) = | ^ | . (41) 



case (cj, 



JV(0 = I C_f <» ) , JVW = ( ^ I ; (42) 



case (d), 



Nl,> = 1 a« 1 1 ' (43) 

Finally let us note that though fermions on two-dimensional sphere without time coor- 
dinate obey to the Dirac equation (|22|) . they are defined by Lagrangian different from dHJ) 
because in the case of two dimensions there are no boost transformations: 

L = y^G~ s (iN + T? 2) e a a D a N + 9 -(<$> 2k N + C {2) N* + h.c.)) , (44) 

where 

D/37 

D A = d A + -f-T[ 2) r (2) - ze(k + kT 3 (2) )A a , 

and e a a is defined in the same way as (jSJ). One can show that this Lagrangian produce the 
same equations of motion as ff22l) . 
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number of solutions 


eigenvalue of 


component 


K < 


-k 


K = —k 


—k < k < k 


K = k 


k > k 


K° 


^0 a 


el 








k + k 


2k 


k + k 


1 


1 


*? 


— K 


-k 










-1 


-1 




k 


— K 


2k 


k — k 






1 


-1 


»/2 












k — k 


-1 


1 



a K a , K° are defined in the next section 

TABLE I. The number of non-zero components in case of different n 
III. INDEX OF THE DIRAC OPERATOR 

In this section we will calculate the index of the operator/} defined in section [TT] by 
Eq.f lTo^ . Indices of Dirac operators without Majorana-like interaction have been considered 
previously (see e.g. 27|, 28] ). The case of the Majorana-like interaction is more involved 
because the Dirac o per ator mixes iV and iV*. For the two-dimensional flat space this problem 



ISO _ T i r 2 

~ lL (D) (D) 



has been solved in 19| . We consider the case of a sphere. 

Let us introduce operators of chirality which distinguish left-handed and right-handed 
modes. In the case g = (no Majorana-like interaction) they are: 

diag (a 3 ... a 3 ), (45) 
K° = T? D) K°. (46) 

These operators anticommute with the operator D^ in . In case g ^ we consider spinors 
(|T3|) having doubled number of components and want to generalize these operators to save 
their anticommutation with D. To achieve this we introduce an analog of K° defined as 

K° 
K° 



K 



(47) 



The analog of K° is 





/ k° 





k = 








[ 


k° 




' k° 





M 











-k° 



in cases (b), (c) 
in cases (a),(d) 



(48) 
(49) 
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Now let us consider an arbitrary Dirak Hamiltonian D (it does not have to have a form 
considered in this paper) and a chirality operator K anticommmuting with D. The fact 
that {D, K} = implies that D transforms left-handed spinors (eigenvectors of K with 
eigenvalue equal to 1 ) into right-handed ones (eigenvectors of K with eigenvalue equal to 
-1 ) and vice versa because 

(1 + K)D = D(1-K), (50) 
(1 - K)D = D(l + K). (51) 

Let us introduce 

Du - D^l, ,52) 

D L = D^JSl. ,53) 

It is easy to check that if D and K are Hermitian, Du and Dl are mutually conjugated: 

D+=D L , (54) 
Dt = D v . (55) 

Therefore, if we have a Dirac Hamiltonian in the space of left-handed fermions, its Hermit 
conjugated operator is the Dirac Hamiltonian in the space of right-handed fermions. 
Now let us define the index as 

index = dim(kerD L D l/ ) - dim(ker D V D L ) = tre"* (z?c,z?L) - tre"* (z?LD[,) . 

This expression can be transformed in the following way: 



index 



( _+n2i=£. t n2 i+k\ 

= tr ( e tu 2 _ e tu 2 ) = 



tr [e- tD21 -^- - e- tD2 ^J = -tr^K). (56) 

Here we made use of a general property of any chirality operator K 2 = 1. 
To find e~ w2 we consider an equation 

^-G(x,y,t) = -D 2 G, G^5(x-y) at t^O (57) 

which is known as the Heat equation; G is its Green function. According to our notations, 
in case g = G is a 2 D I 2 *2 D / 2 matrix and in case g ^ it is 2 D / 2+l *2 D / 2+1 matrix. 
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To continue the calculation of G, we introduce a scalar product in the space where D acts. 
In sec. [Til we defined this space as a set of spinors having a form ffT3]) . Now we consider 
D in the linear span of fundamental solutions of (TT5|) with commutative coefficients. We 
have to do this because if the coefficients were anticommutative, the scalar square of any 
element would be equal to zero. Certainly the change of coefficients from anticommutative 
to commutative does not affect the equations since they are linear, and hence does not affect 
the index. A basis vector in this space is: 



In the case g = the basis element is N®. The scalar product is defined in a usual way 



(58) 



( X m = J dV X + V. (59) 
In case g ^ we define the scalar product as 



(X|*> = J dv( x + *+) ( X j . (60) 



Here dV is an element of volume in S 2 . 

The terms in the expression for the index can be transformed in the following way: 



tre"^ 2 = ^(n | e~ t32 K | n) = 
i 

j dVxdVyY^iN^UxyGij&y^KjkiN^kiy) = 
J i 

J dV x dV y G tJ (x,y)K :jk 5 ik 6(x-y) = 
J dV x G ij (x,x)K ji . 

Here | n) is a basis vector which coordinate representation is (N">) l (x), I numerates the 
basis vectors and i numerates the spinor components of each vector. Also we assumed here 
that 

J2(N^) k (y) + (N^Y(x) = 6{x - y)S ik (61) 
i 
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For g = it can be reached in a straightforward way. Let us consider the case g 7^ 0. If we 
use the representation ffl3|) and note that x = {9 X , (p x } the basis spinors are 





H 


\ ^nmp\ x ) ) 





e p P™(cos s )e im <^ 
e n P™(cos6 r )e- im ^ 



(62) 



Here e p are basis spinors whose p-th component is equal to 1 and the others are equal to 0; 
P™ are Legendre polynomials (In this consideration our basis does not have to be a solution 
of the Dirac equation). 
In this basis, 

£«,„/(y)*UM = m - o.) ( s{ ;< ~ *•> + ) " (63 ) 

nmp ydWx + Vy) 0{f x -lfy) J 

and 

/ 2 D/2+l 2 D/2 \ 

tre- w2 = / dvA J2 G iA x > x ) + ^( G i^^ +i {x,~x) + G 2D/ 2 +i ^{x,-x)) K. (64) 
\ i=i i=i J 

Later (see eq. (176!) ) we see that D 2 is block-diagonal and therefore (see (1C12I) ) G is block- 
diagonal. Consequently, we still can use the expression flBTl) . Finally, 

index = J dV x tr(G(x, x, t)K) (65) 

Now let us find the index in the case g = 0. The expression for D\ in is: 
n 2 _ 

u kin — 

-4 f % + 4s + «* »* - ^ (t ty\ + + c ) , (66) 

/v* \ sin 6/ sin 6/ v ; v ' sin 6/ J 

where 

2 sin 6 1 4 sm & \ v v 

r W^) ( Z |^ + ^))^ + £^ + • (67) 

In derivation of this expression, we took into account that depends only on 9 and does 
not depend on (p. From formula flC12j) we can find the expression for the Green function of 
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the operator (— D 2 ). In case x = y it is equal to 

G(x,x,t) = (I + ~^(e 2 (k 2 + k 2 + 2kKT? D) )A v + 



^-rf^rf^fc + «rf 0) )^cos^) + c 



sin' 6 



Since we know G, we can now easily find the index by making use of Eq. (1651) . The index 
assosiated with the chirality K° is equal to 

index = [ d6dipR 2 sm6tr(G(x,x,t)K) = 2 D/2 [ ^^ekdgA^. (69) 
J Js? 47r 

Since A v is regular in all points except the South Pole, we can rewrite this expression as 

2 D / 2 ek f 

index = / dx a A a . (70) 



4vr j L 

The last integral is taken over a small curve L surrounding the South Pole . We can calculate 
it using the properties ( J2|4l) . So finally the index assosiated with K° is equal to 

2 D ' 2 - l k. (71) 

The index assosiated with K° is calculated in the same way and it is equal to 

2 d ' 2 - 1 k. (72) 

Let us note that by making use of the equations for the vortex, we can transfrom Eq. 
(1691) to express the index through the scalar field. The Lagrangian (CQ) implies the following 
field equations: 

\bd\ ( /~t \ab 



d a (^G s F cd (Gs) ac (Gsr) = (G s ) ab {ie<& + d a <$> - ie<$>d a <$> + + 2e l A a \§\% (73) 



or 



Ae = (Gs) t A(^(G s r(Gsn _ - (74) 



Substituting this expression for A into Eq. (TTUj) . we obtain 



2 D l 2 ik f , n $+d a $ - $<9 a $ 



in 



index = —5—- / dx a ° (75) 
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The index assosiated with K° is obtained from this expression by the change of k to k. 
Now let us consider the theory with g ^ 0. Like the previous case, we find D 2 : 

•>2 _ 



(76) 



^L-^I^I^Cp)^)*^ o 

o {D* kin y - 9 2 m' k (r o {D) rc* iD) T° {D) ( 

In this expression, non-diagonal elements are equal to zero due to Eq. (jXl2j) . In addit ion 
one can notice that T°C(r°)*C* is proportional to a unit matrix. We take into account that 
txK° = trK° = and conclude that this term does not contribute to the index. Thereby 
we have the following expression for the index assosiated with K: 

index = J dV x tr(G{x, x, t)K° + G(x, x, t)*K°) = 2 D/2 k. (77) 

The index assosiated with the chirality K is calculated in the same way. In cases (b), (c) it 
is equal to 

/ dV x tr(G(x, x, t)K° + G(x, x, t)*K°) = 2 d ' 2 k (78) 
and in cases (a), (d) it is equal to 

The result in Eq.f J79|) can be understood in the following way. By looking at the definition 
of K (149]) . we see that if this operator acts to a fundamental solution having a form of (|58|) . 
the obtained vector does not have a form of fl58|) . For this reason we conclude that in these 
cases the zero modes are not chiral. 

Let us emphasize again that we obtained the index in the case g ^ twice as much as 
in the case g = only because we calculated them in different spaces. In the first case the 
space was formed by and in the second case the space was formed by columns (|58l) . If 
one calculates the index in the case g = using the basis f )58|) the obtained answer is the 
same as in case g ^ 0. Actually this calculation can be done just by taking g — in the 
operator D. 



IV. CONCLUSIONS 



We have considered the Dirac equation with Majorana-like interaction on the vortex on 
a sphere. We have found the structure of its solutions which is given by fl39|) . Also we have 
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calculated two indices of this Dirac operator, which correspond to two its chiralities. In 
the case g ^ they are given by (1751 179]) . In the case g = they are given by (I7T| I72|) . 
In principle in high-dimensional space one can consider some other chiralities and calculate 
other indices assosiated with them, but we leave this idea for the further job. 

Now we consider an example of the case d with D = 6. In section [Til the structure of the 
zero modes is written. By looking at the table [TJ one can check explicitly the index for the 
case g = 0. In the case g ^ one can check the index assosiated with K by looking to the 
structure of the solutions described by (1551 -H2D and the list of non-zero components. Also 
one can check that the zero modes are not the eigenvectors of K . 

The last fact is important from the phenomenological point of view. One can consider a 
model where the observable world is M and extra dimensions form S 2 (see ref. 21], 23|). 



In this model K is four- dimensional chirality. So if one assume that the zero modes of N are 
observable particles, they are non-chiral from the four- dimensional point of view, but chiral 
from the extra-dimensional point of view. This is the difference from the particels formed 
by Dirac-like interaction, which are chiral from the four-dimensional point of view (see 21|). 
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Appendix A: Notations 

The manifold considered in this paper is S 2 x M, and M is the Minkovski space. All 
coordinates of this manifold are numerated by capital Latin letters: A, B, . . . .They run 
from to D — 1. 

The coordinates on the sphere are x 1 = 9 and x 2 = (p. 9 is a polar angle and it is counted 
from the North Pole; <p is an axial angle. These coordinates are numerated by small Latin 
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letters: a, b, ... . So the metric tensor of the sphere is 



(GsU = R 2 \ 1 ° o I (Al) 



sm z 9 

Its determinant is designated by G$'- 

G s = det(G s ) ab . 

The coordinates of M are numerated by Greek letters: /x, u, . . . . The zero coordinate x° 
is time. Spatial coordinates are numerated by Latin letters: i, j, ... .The metric tensor for 
the whole manifold is 

G AB = diag( 1 -(GsU -1 1 ) (A2) 



In the case D = 3 gamma-matrices are defined as: 

r (2) = °3, r[ 2) = -ia 2 , r( 2) = io x . (A3) 

The charge conjugation matrix is C = Here a 1 , cr 2 , cr 3 are Pauli matrices. In this paper 
we also consider = iF^F 2 ^ as chirality operator. 

Let us construct gamma-matrices in higher-dimensional spaces. If the number of dimen- 
sions is even and equal to D, we can define gamma-matrices by using (D— 2)-dimensional 
gamma-matrices 7(£>_ 2 ) m ^ e following way: 

T? D) = [ ° r ^~ 2) 1 , A = 0,...,£>-2; (A4) 



D-l 



/ 



r 

L (D) 



1 (D-2) 

-1 

1 



(A5) 



If we want to consider the case of odd number of dimensions which is equal to D + 1 we 
also introduce 

r &) = *(o (A6) 

Also we introduce a projection operator f ^ = diag ^ 1 , — 1 j . One can check that 

rg» = -irg» = H)* +1 r?i» * ... * r^ 1 . (A?) 
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One can notice that our gamma-matrices are imaginary if A = 2, 4, . . . , D(any even 
number more than 0) and real otherwise. In the D-dimensional space we can define charge 
conjugation matrices as 

C(D) = ^Id)^(d)--^{d)^ (A8) 

C (D) = r ?D) r ?D)--- r (D) 2 - ( A9 ) 

In the D + 1-dimensional space only is suitable. 

Defined in this way, C-matrices have the following properties in D > 4: 

rfotfg!) = "^(^(rfz,))*, (Aio) 
rf^S = (-i)^S(r A r- (ah) 

Matrix Cp) written without upper index denotes any of these charge conjugation matri- 
ces. 

By making use of the last formulae, one can show that 

r (D) r fD) r ?D)C(D) + r°( D) C( D) (r°( D) )*(rf D) )* = o. (Ai2) 

This relation is useful in derivation of many of the expressions mentioned in this paper. 
One can show that 

(^^^(-l)^- 1 ^)^), (A13) 

(O^h-i)^- 1 ^)' 

(C$r = (-l)- ( - +1) ^), (A14) 
(<$))'<$) = (A15) 
(<$))'<$) = (-1) ( - +1)( -- 1} . (A16) 
Also in this paper we use the properties: 

{Tf D) ,T° {D) C (D) } = if D/2 mod 4 = 1 or 3 

[Tf D) ,T° (D) C (D) ] = if D/2 mod 4 = 2 (A17) 
[Tf D) ,C (D) ] = if D/2 mod 4 = 1 or 3 

{Tf D) ,C (D) } = ifD/2mod4 = 2 (A18) 
The structure of C-matrices in the cases used in this paper is the following: 
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(a) D/2 mod 4 = 1, 



C# n. o 

-c, 



C S = M „ ^> (A19) 



(D-2) 



(b) D/2 mod 4 = 2, 



^. I J! °^- 2) I (A 20) 



'(©-2) 



(c) D/2 mod 4 = 2, 



°(X>-2) U 



(d) D/2 mod 4 = 3, 



= i c< r H ° (A22) 



Appendix B: The solution of the equations for zero modes 



Here we explore the fundamental solutions of the equations (1361 - 1371) . In the case g = 
we find them explicitly and in the case g ^ we find their asymptotics. We assume k > 0. 
Let us write these equations in the component form: 



(de + te(n + k)A 9 )e n - igR^vh+i-n = 0, 

2 smv 

cot i 

(d e + —)£ + — ^(^ - <k + fc)^)& + ^$ 2fe ^ +1 _ n = 0, 

(d e + ^)fo]L ) + -r^a " - k)AM +1 _ n ) + igR^f^l = 0, 
2 siii0 

Pot r/ ? 

(do + -^-)(v 2 N + x- n ) ~ ~ «e(* - k)A v )(rf N+1 _ n ) - igRm 2k e n = 0; (Bl) 
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pA 
sin 9 



sin 9 



To get rid of the dependence on ip we use the following ansatz: 

^1 = u L exp + + jf 

^ = Mc/ exp N(Z + \)^ + J q 
Vn+i-u = Vl exp ^-i(2fc + Z + ^ + jf ^^(~« + , 
^ +1 _ n = ^exp(^-z(2fc-Z--)^ + y d0_|( K -A:)J. (B2) 

Here «£, ujj, vl, vjj are new variables depending only on 9, and Z is an integer number which 
further will numerate the solutions. Let us note that in this ansatz, we assumed that the 
spinor components obey antiperiodical boundary conditions. They are chosen because if we 
move on a closed curve from (p = to if = 2tt we rotate our velbain (JSJ) by 2ir. This results 
in the change of the sign of a fermion. In terms of the new variables the equations read 

1 + 1/2 COt 2fc o r e K A v de 

Oqu l ^rU L + —u L - gRf e J iv L = 0, 

sin 9 2 

1 + 1/2 COt ok 2 f eKA f dB 

oeuu ^-uu + —u v + gRf e i iv v = 0, 

sin 2 

2fc + Z + 1/2 cot nr2t „ f f^V!. 

<9^l + . „ v L + —v L + gRf 2k e 2 J **eiu L = Q, 

sin 6* 2 

_ 2k -I -1/2 cot nf2 (. orf^V!. rt 

. 5 v v + —v v - gRf 2k e 2 J lUv = 0. B3 

sm 2 

Then we make another change: 

(sin 0/2)' 



«l = Vl 



v l = z L - 



uu = Vu 



COS9/2Y+ 1 ' 
s\n9/2Y 2k - 1 - 1 
(cos9/2)- 2k - 1 
(sm9/2) 1 



{cos9/2) l + v 
{sin9/2) 2k - 1 - 1 
(cos#/2) 

to we obtain the following equations: 



<% L - ^ J R/ 2fc e- 2 ^ f ^(tan^/2)- 2fe - 2i - 1 ^ = 0, (B5) 
d z L + igRf k e 2 S^(tan9/2) 2k+2l ^y L = 0, 
deyu + igRf k e 2r -^(tan9/2) 2k - 2l - 1 z u = 0, 

- igRf 2k e- 2 f^ l (t8,n9/2)- 2k+2l+1 y u = 0. (B6) 
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Now let us consider separately the case g = 0. Here we conclude that all variables in the 
equations flB6j) are constants: 



yL = y° L = const, 
yu = Vu = const, 
Zl = z° L = const, 



zu = Zjj = const. (B7) 
Substituting them into the expressions (1B2I . IB4j) we find the solutions for the initial variables: 



„ 2 (sin 0/2)' / 1, f e 7n eA v . 



- ^ ty ^ (~ l{2k - 1 - & + l° M ^ (K - k) ) ■ (B8) 

These expressions may have singularities at 9 — > or 9 — > ir. To investigate them we write 
the asymptotics of Eqs. ( 1B8I) which are the following (£ = it — 9): 

• 9^0: 

e n ~yUe/2) l ewUl + l)v 



9 ->• 



7r: 



^ + i-n ~ z° L (9/2)- 2k - 1 - 1 exp f-i(2* + i + \)^\ , 
»&H-» ~ ^/a) 2 *-'" 1 exp f-i(2* - J - i)^ . 

V 2 N+ i- n ~ z° u C K - k+l exp (-i(2k - I - i)^ . 



(B9) 
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These expressions are regular if 9 and ( have non-negative powers. From this condition 
we find the constraints for 

£n : < I < ~k- k- 1, 
&:0<l<k + K-l, 
Vn+i-u :-K-k<l< -2k - 1, 
VN+l-n :K + k<l<2k-l. 

Finally we can write the structure of the solutions depending on k and k (recall that 
k > 0): 

K + k < — 1: £^ has \k + k\ modes and rf N+l _ n has k — k modes. 
K + k = 0: Only rf N+l _ n has non-zero modes. Their number is 2k. 
\<K + k<2k — 1: has k + k modes and rf N+l _ n has k — k modes. 
k + k = 2k: has 2k non-zero modes. 

K + k > 2k + 1: nas K + ^ modes and rj l N+1 _ n has k — k modes. 

Now let us consider the case g ^ 0. A general analytical solution to the equations ( 1B6|) is 
unknown. In this paper we find the asymptotics of the solutions in cases 9 — > and 9 — > ir. 
To do it we derive the second-order equations containing only y L , yu from Eqs. (IB6|) : 



2 fl + 2k + 2l d e f 2enA\ ( ^2Mk n 

+ r-r 2 A;— + d e y L - {gRyp k y L = 0, 

\ sine/ / sinS / 

2 f l-2k + 2l dgf 2enA\ 2 4fc 

\ sm# / sm9 J 

Each of these equations has two linearly independent solutions. After substituting the 
asymptotics of the vortex fields, eq. (OH), we can find the behavior of the unknown functions: 

• 9 0: 

y L ~ 2 z c ( 1 + , , g , °f , -fl 4fc+2 + 2%g9- 2 \ 

y °V 2(4fc + 2/ + 2)(2A; + l) J oy 

(„2 T32(~iAk \ 
2(2Z + 2)(2fc + l) / y 
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• 9 -> 7r: 



yL 2 C - ^"4(fc + K + /)J +2 ^ C 

W7 ~ 2^ (l + ff^ ) + r^MC 2 '-^^. (Bll) 

Here ao, &o> c o, cfo, a 7r, &tt, cv, are arbitrary constants and factors 2 l , 2~ l ~ 1 and g are intro- 
duced here for convenience. 

If we substitute the last expressions into Eqs. ( \B6\ IB2I . IB4j) . we find the behavior of our 
initial functions: 



il - (c e l + d og e- 1 ) e-^ +1 /% 

, / -ZC ^i?C 2 ^6' 2A:+ ' +1 zdp(2j) 2fc _;_A i( 2 fc+Z + l/2)^ 



4A; + 2/ + 2 ' i?C 2 £ 
igRCl k F a Q 6 2k+l+l ib (4Jfe - 2/) 



2 / ^Jto 0F u, iy mtq^ - *t) Q2k -l-l -i(2k-l-l/2)<p 

VN+1 ~ n ~ V 2TT2 + Scg J 

tt: 



( -igRCf F c^- k ~ K - 1 , ^(2/ + 2fc + 2/, + 2)C ;+fc+K ^ _ w+1/2)v 



F 



C 



To select regular solutions we use the so-called different signs theorem (see also [21] . 



201]) . It states that if y^u and dei/L,u have the same sign at one interior point #0 of the 
interval (0,7r), then they have the same sign at all points of the interval (0q,tt). Let us 
prove it. If a function and its derivative have the same sign which is positive(negative) the 
function cannot change it before its derivative does it. The derivative can change its sign 
in a point of a local maximum (minimum). From Eqs. (1B10j) we see that if dgyi,u — and 



dgf/f, A v do not have any singularities (this is true for the vortex fields), then y u>L and 
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9gVu,L have the same sign. Thus we got the contradictory statement that in the point of a 
local maximum (minimun) the second derivative is positive (negative). 

Now let us consider the solutions £^ and ^ +1 _„. One can see that their linearly inde- 
pendent components have opposite powers of ( in case 9 — > n so only one of them should 
be non-zero if we want to get a regular expression. So there are two possible cases: 

— k — k, — I — 1 > : c n ^ 0,d n = 0, 

k + n + l>0: ^ = 0,4^0. (B13) 

Since the asymptotics have the power-like form, we conclude that in both cases if £ increases, 
Vl(C) also increases. If we want £^ to be regular, only one of the coefficients c ,d should be 
non-zero. If we suppose that do ^ 0, then from the expression for r/^ +1 _ n we conclude that 
/ < —2k — 1 and therefore Vl(9) increases. This contradicts to the different signs theorem. 
If we suppose that Co ^ 0, from the expression for ^ we conclude that I > and so yi(9) 
also increases if 9 increases. This possibility also contradicts to the different signs theorem. 

The case of and rj 2 N+1 _ n is treated in the same way. From the asymptotics of 9 — > 7t 
we see that there are two cases: 

k + K-l-l>0: 0^^0,6^ = 0, 

l-k-K>0: a n = 0,b 7r ^0. (B14) 

In these cases (as opposed to the previous ones) we concede that ao,&o 7^ together. 
From the condition that rf N+l _ n cannot have negative powers of 9, £, we conclude that 
< I < 2k — 1. Also, if we force these coefficients to have different signs, we can satisfy 
the different signs theorem. Let us consider the case when only one of the coefficients do, bo 
is non-zero. For ao ^ we have I > and so yu{9) does not satisfy the theorem; and for 
bo 7^ we have I < 2k — 1 and yu(9) again does not satisfy the different signs theorem. 

Thus we have shown that and rf N+l _ n have regular solutions but £^ and ^ +1 _ n do 
not have them. 

One can note that the solutions in the cases g = and g ^ have different structure. 
Nevertheless, if we take the expressions (1B12j) and find their limit in case g — » 0, we obtain 
the expressions having the same form as Eqs. flB9l) . 
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Appendix C: The Green function of the Heat equation 

Let us consider a general Heat equation 

j f G(x, y, t) = {a^d a d p + b a d a + c)G(x, y, t) (CI) 

supplemented by the condition 

G(x, y, t) -> 6(x - y) at t 0. (C2) 

Here d a = and a, b, c depend on x a . We can solve this equation in the case x ~ y and 
t — > with the following ansatz: 

G(x,y,t) = e s ^ t Y,M*,y)t k . (C3) 
fc=-i 

By substituting it into Eq. fIClj) and equating the terms having equal powers of t, we get 
the following equations: 

S = a af3 d a Sd^S, (C4) 

A_i = -a^dadpSA-i - 2a a ^d a Sd /3 A^ 1 - b a d a SA_ u (C5) 
= a aP d a d p 3Ao + 2a a/3 <9 a S<9 /3 ,4 + 6 a 9 Q SA + 

{a aP d a d p + b a d a + c)A- V (C6) 

Now we take into account that x ~ ?/ and assume a, 6, c(x) ~ a, 6, c(y).The equation (IC4|) 
has the solution 

S = --a afl Ax a AxP, (C7) 



where 



a Q/3 = (a M ) , 



Ax a =x a - y a . 
Using the expression for S we transform the Eq. (1C5|) to 



It's solution is 



Ax a d a A^ + baa ^ Axl A _ x = o. (C8) 



A_! = e ^ (C9) 
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By substituting S and A_i into Eq. (1C6j) . we get 

= -Aq - Ax^^Ao ^ A + ( ^ + c)A_i, (CIO) 



Its general solution is 



• b a Wa a a const , , 

A)= r^ + c +^^^ ^-i ( cn ) 



2 

a V 



To avoid a singularity in the last expression we take const = 0. 
Finally after substituting (1U7|) . (jU9]) . (iCTTj) into fl£3J) we have 

G = - f i - + c) e -h^^-^¥^ (C12) 



4tt \t 4 

The constant j- is introduced here to satisfy the condition flC2j) . 



[1] N. Arkani-Hamed, S. Dimopoulos, G. R. Dvali and J. March-Russell, Phys. Rev. D 65 (2002) 

024032 |arXiv:hep-ph/9811448| . 
[2] V. A. Rubakov, Phys. Usp. 44 (2001) 871 [Usp. Fiz. Nauk 171 (2001) 913] 

[arXiv:hep-ph/0104152] . 
[3] V. A. Rubakov, Phys. Usp. 46 (2003) 211 [Usp. Fiz. Nauk 173 (2003) 219]. 
[4] Kaluza T Sitzungsber. Preuss. Akad. Wiss. Berlin, Math.-Phys. Kl. (1) 966 (1921) 
[5] Klein O Z. Phys. 37 895 (1926) 

[6] V. A. Rubakov and M. E. Shaposhnikov, Phys. Lett. B 125 (1983) 136. 



[7] M. Visser, Phys. Lett. B159 (1985) 22. |hep-th/9910093|. 

[8] R. Jackiw and C. Rebbi, Phys. Rev. D 13 (1976) 3398. 

[9] R. Jackiw and P. Rossi, Nucl. Phys. B 190 (1981) 681. 
[10] Klaus von Klitzing, The quantized Hall effect, Nobel lecture, 1985 
[11] L. Fu and C. L. Kane, Phys. Rev. Lett. 100 (2008) 096407. 

[12] Y. Nishida, L. Santos and C. Chamon, Phys. Rev. B 82 (2010) 144513 [ arXiv:1007.220Tl 
[cond-mat.supr-con]] . 

[13] J. K. Pachos , M. Stone and K. Temme, Phys. Rev. Lett. 100, 156806 (2008), 

larXiv:0710.0858k r3 
[14] A. Kitaev, 2001 Phys.-Usp. 44 131, |arXiv:cond-mat/00i0440l 

29 



J. K. Pachos, Cont. Phys. 50, pp. 375-389 f2009^. larXTv:0812.1116l 

A. Abrikosov, ZhETF (USSR) 32 (1957) 1442 [JETP (Sov. Phys.) 5 (1957) 1174], 

H. B. Nielsen and P. Olesen, Nucl. Phys. B 61 (1973) 45. 

E. Witten, Nucl. Phys. B 249 (1985) 557. 

E. J. Weinberg, Phys. Rev. D 24 (1981) 2669. 

M. V. Libanov and S. V. Troitsky, Nucl. Phys. B 599 (2001) 319 |arXiv:hep-ph/0011095| . 



J. M. Frere, M. V. Libanov, E. Y. Nugaev and S. V. Troitsky, JHEP 0306 (2003) 009 
[arXiv:hep-ph/0304117] . 

M. V. Libanov and E. Y. Nugaev, Surveys High Energ. Phys. 17 (2002) 165. 
J. M. Frere, M. Libanov and F. S. Ling, JHEP 1009 (2010) 081 |arXiv:1006.5T96 [hep-ph]]. 
M. Hirsch, S. Morisi, J. W. F. Valle, Phys. Rev. D78 (2008) 093007. |arXiv:0804.152T1 [hep- 
ph]]. 

A. Roy and M. Stone, J. Phys. A 43 (2010) 015203 [arXiv: 0909 .15691 [cond-mat.other]]. 
T. T. Wu and C. N. Yang, Phys. Rev. D 12 (1975) 3845. 



D. V. Vassilevich, Phys. Rept. 388 (2003) 279 [arXiv:hep-th/ 0306138 1. 

A. S. Schwarz, "Quantum field theory and topology," Berlin, Germany: Springer (1993). 

R. Jeannerot and M. Postma, JHEP 0412 (2004) 032 |arXiv:hep-ph/0411259| . 

M. B. Green, J. H. Schwartz and E. Witten, "Superstring theory," Cambridge, UK: Univ. Pr. 

(1987) (Cambridge Monographs On Mathematical Physics ). 

V. A. Rubakov, "Classical theory of gauge fields", Princeton University Press, (2002). 
K. V. Stepanyantz, "Classical field theory", Moscow, Fizmatlit, (2009), (in Russian) 



30 



